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We found the contribution to the vacuum expectation value of the energy-momentum tensor of 
a massive Dirac field due to the conical geometry of the cosmic string space-time. The heat kernel 
and heat kernel expansion for the squared Dirac operator in this background are also considered 
and the first three coefficients were found in an explicity form. 



I. INTRODUCTION 



The study of quantum fields in the space-time of a static cylindrically symmetric cosmic string has been done in 
different situations. Examples of these studies are the computation of the non-vanishing contribution to the vacuum 
expectation value of the ener gy- momentum tensor of quantum fields, as for example, scalar, spinor and vector fields 
[2I H. Ull 13 H7L J23I I25L 13^. Jill I44I liBl Eil . Quantum fields at non-zero temperature were also considered in the 
literature [3, flCM . Ha . l24l l35l 1361 l37l l38l |43 | . It was noted that the energy-momentum tensor has a non-i nteg rable 
singularity at the origin and possesses the structure similar to the energy- momentum tensor in the wedge |l3j . On 
the other hand there is a great interest connected with properties of the heat kernel of the Laplacian like operator 
and its expansion in this background. The point to be considered is that in this space-time there is a singular surface 
with codimension two which originates some problems [33|. As a consequence, it is impossible to apply standard 
expressions for the heat kernel coefficients as discussed in |46j and therefore we have to modify them. Nevertheless, 
there is a way to calculate the heat kernel coefficients by considering generalized cone in dimensions greater then four 
and to treat the heat kernel coefficients in the standard way using dimensional reduction. General considerations 
concerning the first heat kernel coefficient and up to systems with spin-2 were considered in Ref. |'2C| . 

Cosmic string j^, 01 is an exotic topological defect which may have been formed at phase transitions in the 
very early history of the Universe. Up to the moment no direct observational evidence of their existence has been 
found (see, nevertheless H3), but the richness of the new ideas they brought along to general relativity seems to 
justify the interest in the study of these structures. 

The gravitational field of a straight infinitely long cosmic string is quite remarkable; a particle placed at rest around 
a straight, infinite, static string will not be attracted to it; there is no local gravity. The space-time around a cosmic 
string is locally flat but not globally. The external gravitational field due to a cosmic string may be approximately 
described by a commonly called conical geometry. Due to this conical geometry a cosmic string can induce several 
effects like, for example, gravitational lens 0,E3, P air production |2fj, electrostatic self- force |3_ on an electric 
charge at rest, bremsstrahlung process and the so-called gravitational Aharonov-Bohm effect |lr]|. 

In the present paper we obtain an explicit form of the expression for vacuum expectation value of the energy- 
momentum tensor of a massive Dirac field in the cosmic string space-time, showing up that the conical structure of 
this space-time induces a non- vanishing value for this quantity. The massless case was already considered by Frolov 
and Serebriany |T^ . The energy-momentum tensor is traceless and it has the same non-integrable singularity at 
the string as in the massless case. We also reobtain the heat kernel and heat kernel expansion in terms of contour 
integrals. We show that in the coincidence limit, the Green function of squared Dirac operator is defined by the first 
heat kernel coefficient and the energy-momentum tensor is defined by the second one. The vacuum expectation value 
of the energy- momentum tensor of a massive scalar field was considered in Refs. p3ll2^|. 

The paper is organized as follows. In Sec. UTI we find the Green function of the squared Dirac operator and the 
corresponding heat kernel. In Sec. IHII we obtain the expansion of the heat kernel and found in manifest form the heat 
kernel coefficients which are expressed in terms of the delta function and its derivatives. The vacuum expectation 
value for the Dirac field is obtained in Sec. IIVI We conclude with the Sec[V]by summarizing the results obtained 
presenting some remarks. In this paper we adopt the units where h = c = G = 1. 
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II. GREEN FUNCTION AND HEAT KERNEL 



To start with let us consider the massive Dirac equation in the Euclidean sector, r = it, in the space-time generated 
by a cosmic string. The space-time is described by the following line element |47| 

P 2 

ds 2 = dr 2 + dp 2 + hrdip 2 + dz 2 
v l 

where r, z £ K, p € M + and <p £ [0, 2ir]. The manifold associated with a cosmic string, M. = M 2 ®C 2 , is a direct 
product of two the dimensional space K 2 and two dimensional conical space C 2 . The scalar curvature has to be 
understood as a distribution 431 



K = 4tt 



v- \^ 2 \x) 



(1) 



The general considerations about the properties of this space-time may be found in Q while the Green function of the 
Dirac equation, for a massive field in this background was found in Ref. |37j . We also obtain the Green function in 
this background, but in different form compared with the one obtained in [37| which is more suitable for our proposal. 

According the general considerations concerning the Dirac equation in arbitrary dimensions 0, 122I l32l | , the Eu- 
clidean spinor Green function obeys the following equation 

P m S(x;x') = -h y — 

V9 

where TJ) m — 7 M V M + m is the Dirac operator and V M = + T fl . The Fock-Ivanenko spinor connection 

i M=2 CT e (a) V M e (b)^ 

is expressed in terms of the vierbein e(M„ and the matrix a ah — j[y^ a \ 7^]- Gamma matrices satisfy the Clifford 
commutation relation { r f^ a \ 7^} = 25 ab Ii and the generalized Dirac matrices, 7 M (x), are defined by the relation 
7 M (x) = G^ a ^ a \ In what follows we will adopt the standard representation 



v(°) 



1 

-1 



a k 
-<r k 



where a k are the Pauli matrices. As to the vierbein, we will consider the one used in |l2j . which are given by the 
following expressions 



(\ 



-(a) 







cos if — j sin ip 
sin ip +-cos(p 



\0 







y 



where rows correspond to the number (a) of the vector ef Thus, we obtain the results 



7 (0), 7 



p _ 



cos <p7 



(i) 



- sin ipj 



(a)" 

(2) ~V - 



7^ = sinyyy^ H — cos ipj^' , 7 

P P 



(2) 



and 



r \~v>~Q)$ =i- -s„s 



2v 



2v 



where 



i-y 7 . Therefore, the spinor Green function S obeys the following equation 



v- 1 



§*( x _ rf) 



Let us define the Green function £/ of the squared Dirac operator by the relation 

S(x;x')=JZ>_ m g(x;x'). 



(2) 
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where 



It obeys the following equation 

ft G{x;x ) = -I 4 — , 

p 2 = P m p_ m = [sT^ + l -d p t=£? + i^^-E 3 d v -m 2 - Jw) . (3) 

This equation has no the form of the equation for the scalar Green function. Its general form was presented in |46| . 
In this paper we will omit the term with the singular scalar curvature. The contribution arising from it has been 
discussed in the literature [3. 13, IM l29j . 

Due to the fact that S3 is diagonal, the Green function Q is diagonal, too. In order to obtain the Green function 
in explicitly form let us find the full set of bispinors which satisfies the equation 

p 2 <j>{x) = -A 2 0(x). (4) 

We shall numerate the eigenfunctions by eigenvalues of operators p T , p z and the projection of the full momentum J3, 

as 

p T <f> = —id T <t> = E(f>, 
Pz<P = -idg<f> = p s <f>, 

^)<f> = {-id, + \ 



■h4> = (p v + ^-S 3 )0 = {-id v + -S 3 )0 = j4>, 



and define the number p±^ by relation A 2 = p 2 ± + E 2 + p 2 + m 2 . All these operators commute with each other and 
with Dirac and squared Dirac operators. The above set of operators define the eigenfunction up to four arbitrary 
constants. Let us numerate these independent solutions by eigenvalues of matrices £3 with eigenvalue a = ±1 and 
75 = 7^°) 7W <y( 2 ) <y( 3 ) with eigenvalue b = ±1. Both these operators commute with the above set of operators. Wc 
shall numerate eigenfunctions by numbers q = E,p 3 ,p±,l where E,p 3 € K,pi S K + , I = 0, ±1, . . . (J = I + |) and by 
(a, b) = (±1,±1). Therefore we have the following four independent solutions 



</ 

-,±) 



47r 3/2 
47r 3/2 



where 



w (+,±) 

are the eigenfunctions of S3 and 75. Here J M (z) is the Bessel function with indices (3^ = \uj =F i|- This full set obeys 
the following relations of completeness and orthogonality 















loi 


« w (-,±) = 











Y,jJJ dEdp^dp^^{x')^ b \x) = 5{A){X j_ X '\ 



4> ( ^ b,) \x')<t>^ b \x)^d 4 x = 5 a>a ,6 b , bl 5i, v 5(E-E')5( P3 -p' 3 ) 



P-L 



Now we are ready to obtain the Green function in manifest form. To do this, let us represent the Green function in 
the following form 

- EE/// (5) 

a.b I - 1 - ° 

v r r r p iEAr-\-ip 3 Az 

irfL dEdp 3P± dp A 



(2n)*^JJJ ^^pl + E^+pl 
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x diag(J /3 _(r)J /3 _(r')e jiA ^ J 0+ (r) J p+ (r')e^ l+1 ^, J Ar)J p _(r')e aA ^ J l3+ (r)J f}+ (r')e^+ 1 ^), 

where r = p±p, r' = p±p'. By changing the variable I — » / — 1 in the second and forth terms, we can write the Green 
function in more simple form, as 



Q{x;x') 



iEAT+ip 3 Az+UA(p 



V \ III g.J-^^J. I 



' (27 

x diag(J ?+ (r)J 5+ (r'), J e _(r)J s _(r'), J e+ (r)J 6+ (r'), J^_(r)J e _(r')), 



where f± = \vl±^\. 

It is easy to find from this expression the heat kernel of the squared Dirac operator in Eq. J3J. It is given by 



K.{x\ x'\s) 



(2tt) 3 



Y, 1 1 1 dEd P3P± dp ± e iEAT+t P 3Az+aA ^- s ^+ E2+ ri +m2 " 



x diag(J c+ (r)J c+ (r'), J e _ (r) Jf_ (/), J S+ (r) Jf + (/), (r)J^ (/)) 
= /C^(Ar,Az| s )/Cl 2 )(A^,p,p'| s ), 



(6) 



where Km (Ar, Az|s) is the heat kernel of the operator L m = d 2 2 + d 2 2 — m 2 and K,)f (Aip, p, p'\s) is the heat kernel 
of the operator 



"(2), 



L v = d 2 p2 + -d p + V —dl, - 



(v- l) 2 t .v{v- 1) 
4p 2 



■' 2 _S 3<V 



Integrating over E,p 3 ,pj_, we get 
*£>(Ar,A*|.) = 



2s 



2s 



2s 



2s 



iE^-— <imw£>. w|i w£>, tiji) 



2s 



2s 



2s 



2s 



(7a) 



(7b) 



It is worth noting that this expression leads to a heat kernel which is a diagonal matrix. The components of this matrix 
are defined in terms of the first component ICi t i which we denote by K by relations: £3,3 = K and /C2,2 = ^4,4 = K*. 
Indeed, because £± — ► £^ with changing I — > I then 



1 1 V 1 J 



and the heat kernel has the following structure 

K, {2) =diag(K,K*,K,K*). 
Using Eqs. J2J) and l@ we obtain the spinor Green function in manifest form, which can be written as 



ds 



,2 ,-T 2 + Az 2 + p 2 +p /2 



■(_ m5l __I 52 +^3-^) 



where 



Si = diag(l _,I 0+ e iA ^,I _,I 0+ e iA ^), S 2 = diag (l p _ , Ip + e iA *, -I _ , -Ip + e iA *) 
( I _ \ 



S3 



5 4 



-I 0+ e iAt ? 



y I 0+ e iA f / 

/ 



_ e -V( p ' J/3 __ pJ/3+ ) 

\-e^{p'I 0+ ~pI p _) 



e-*{pfIf l _-pIf l+ )\ 

e^(j/l fi+ -plp_) 



y 



(7c) 



(8) 



(9) 
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and the argument of the Bessel functions of second kind Ip ± is z = To obtain previous expression we used the 
relations : 

where e; = sgn(Z). 

Now, let us consider some limiting cases. In the massless case, to — 0, the Green function can be found in a closed 
form and is given explicitly by 

„, „ v 1 1 sinh(^77)/ 4 + sinh(^77)e- lA ^ E3 

S(x;x ) = ^-2— -T-r r 7; • 

»7r^ pp' smh r/ cosh ^77 — cos Aip 

We note that this Green function does not coincide with that obtained in [T7L |37| and it does not obey the relation 

G(p + 2ir) = -g(ip). 

Instead of this, it obeys the relation above with minus sign changed toplus sign. The reason for this fact is connected 
with the use of another vierbein as compared with the one used in [17j . 

In order to reobtain the result showed in |l7j , we have to apply a rotation of the vierbein we have chosen by an angle 
ip, in the plane (p, ip). It is well-known that the bispinor will be multiplied by a factor e^ vE3 . The Green function, as 
a bispinor at point x and point x' , is multiplied by factor e5 A( ^ S3 . Thus, multiplying by this factor we arrive at the 
formula 

, . v 1 1 ei A ^ 3 sinner?) + e ~^ A ^ 3 sinhf^ir?) 

G(x]x) = ^j — ^-r \ r 

57r z pp' smh 77 cosh vr\ — cos A if 

which coincides with the results obtained in 0,H3|, if we take into account that ^£3 = -yW'-yt 2 ). 
In the case of zero angle deficit, v — 1, we have £± = |'| and thus 

where R — (At 2 + Az 2 + p 2 + p' 2 — 2pp' cos Alp) 1 / 2 is the distance between points x and x' in the Euclidean 
4-dimensional space in cylindrical coordinates. The Green function in this case turns into 



f°° 1 TO 

G(x\x')= \ K(x;x'\s)ds = —?—K\(mR)h, 
Jo 47r R 



as it should be. Here is the Bessel function of second kind. 



III. HEAT KERNEL EXPANSION 



Let us investigate the asymptotic expansion of the heat kernel of the squared Dirac operator which is given by Eqs. 
(I7bl) and (|7c|) . It is impossible to take an expansion of the Bessel function because the asymptotic expansion of the 
Bessel function depends on the ratio argument and indices. For this reason we use the integral representation for the 
Bessel function, as follows 

IJz) = — [ e i{w ~^ ] ^ +zsinw dw. 

The contour T lies in the half-strip fcg [— |, 5z > 0. It goes from — f + ioo to ^ + ioo. 
Let us consider, in general, the series 

1 

where |a| < 1. For our cases we have a = and |a| < \. Then we use the representation above for the Bessel 

function and change the integration and summation and shift the variable w to z = w — ^. In this situation, we 
obtain the following expression for the function /(a, r) 
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where e a = sgn(a). The contour V is the contour F shifted by tt/2 to the left. It lies in the positive part of the 
strip z € [— 7r, 7r]. Then we divide the integral in two parts according with integrand and in the second part, we do 
the change of variable z — > — z. The integrand takes the form as the first integrand but the contour is the central 
symmetry of V . Therefore we may recombinc both integrals into a single integral and arrive at the following formula 
for function f(a,r) 



f(<*,r) 



-iz\ua\+r cos z 



2tt .L e -i(.xv-e a Aip) - i 



dz. 



The contour 7 has two branches. First one is contour V and the second is the central symmetry of V . Taking into 
account this formula we arrive at the following expressions for heat kernel component 



K 



2, 11 



4-7TS 



-/( 



'-1 w\ 

2v '2s' 



8n 2 s 



-i(zv—A<p) "Y 



■dz. 



The complex conjugate K* differ from K by the sign of Aip only. 

Now, let us return to the function f(a,r) and modify the contour 7. It is easy to see that the zeros of the 
denominator are in the points z n — e a ^r + Because of \Aip\ < tt and v > 1, then we have that|Ay>/z/| < tt. 
Therefore, the first zero of the denominator, zq — c a — — , belongs to the interval of integration [— tt,tt]. Let us extract 
this first pole in manifest form. With this we add the integrals over two lines around point z — e a Aip/v in each of 
which we integrate in opposite directions. The contribution from these additional integrals are zero. Then we divide 
the contour 7 with additional lines in two parts and represent the integral as an integral over contour x an d over a 
circle around point e a —jf-. The contour \ has two branches. The first one goes from — tt + ioo to — tt — ioo, very close 
to the point e a ^- from the left side. The second part of the contour goes from tt + ioo to tt — zoo, very close to the 



point e a -^r from the right side. The second integral is minus the residue at point e, 
our function /(a, r) in the form below 



Am 



Therefore we may represent 



A 1/ T &~ iz\va\+r cos z 

J V ' > ^ 2tt J x e -<zu-e a A v ) _ l 



dz. 



Taking into account this representation we arrive at the following expression for component K 



K 



Atts 



y p - — ilzis—Atp) '^ 2w 1 +^-(cos2-cos -77-) 

1 + 77- / -f 7T~\ dz 

2tt / v e -*(zv-A v ) _ 1 



where d 2 = p 2 + p' 2 — 2 pp' cos The expression for the complex conjugate K* component differs from the above 
by the sign of Aip. 

Therefore, in the coincidence limit p' = p and ip' = ip we obtain the following expression for the heat kernel 



JC^(x-x\s) 



Atts 



2tt j % 



-dz 



(10) 



where the contour \ intersects the real axis close to the origin. In order to find the expansion of the heat kernel over 
s we use the formula presented in , which is given by 



^ n\{Ao) n ' 



n=0 



where 5^ 2 \r) — ~5(p)8(ip) and A( 2 ) = d 2 2 + ^d p + tt^2 ■ The integral may be calculated by the method of residue 
at the zero point. Thus, we obtain 



4 2) (*;z| S ) = ^ 



n=l 



where 



a n (x;x)=h (n ^ n _ 1 A^(r) 
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and 



R n = Resg 



{e~ izv - 1) sin 



2n+2 



Therefore, the first three coefficients are given by 

Rq - 



n{v 2 -l) 7^ 2 + 17 D 31 f 4 + 178^ 2 + 367 

i -"-I = ™ -fKh Rl = r.,-.^ #0. 



6^ 60 2520 

The expansion of the four dimensional heat kernel is expressed in terms of the same coefficients, as 



JC(x; x'\s) 



1 



(4tts) 2 



h + ^2 a n (x; x)s n 



From the above expression we observe that the first coefficient Rq for the spinor field is minus 1/2 of the same 
coefficient corresponding to the scalar field. Therefore, we have 



tra 



(1/2) 



2 1 



-2a 



(o) 



and thus in accordance with fl9t |20L 1281 . The other coefficients have different structure. 



IV. VACUUM EXPECTATION VALUES 



First of all let us find in closed form the coincidence limit of the renormalizcd Green function of the squared Dirac 
operator. To do the process of renormalization we subtract the same function with v = 1. Because the second integral 
in Eq. (|lf)|> is zero for v = 1, we have the following expression in the coincidence limit 

r°° rfg r g — iz iL ^ L - £ T- sin 2 m 2 s 

32tt 3 In s 2 / ~ < ; " 1 



S ren ^*) = ^=sl -si dzh. (11) 



x 



Firstly, let us consider the massless case. We may integrate easily with respect to s and obtain 



32tt 3 p 2 J x { e - lzv - 1) sin 2 § 

Then we close the contour to infinities and calculate this expression with the residue at the point zero. Thus, we get 
the result 

r ren, \ vRq V — \ 

lbir A p z 9bn /! p z 

In the massive case, we deform the contour in two straight lines over z = ±7r. We have to take into account the 
residues at points z n = 2nn/v, where |n| < v/2. For v < 2 there is no zeros in the denominator and no residues. By 
changing the variable z = ±n + iy we obtain, in general, the result 

e _^_Z si „ 2f np -oo p2 sinhfsinhf 

: dz = 4cos — / e ~~ ccteh 2 - - — dy. 

x e ,z "-l 2 J cosh vy — cos nv 

Taking into account this formula we arrive at the following expression 

„».«./ x m v-?/ x„ / tt^s mi/ cos ^ Z" 00 Xi(2mpcosh f ) sinh f sinh ^ 

£T« x;x = — — ^ -l "tan— Ktfmpsm— )+ - 3 2 / , , 2 r 2 dyh. 

P ti—i v v P J° cosh i| cosn z;?/ — cos 

In the case f < 2, only the last term survives. The appearance of additional terms, when v > 2, is related with the 
fact that additional closed geodesies in C 2 , appears in this case [3(j. For rap ^> 1 the above expression exponentially 
falls down, according to 

qr^(x-x)- vRo m2e ~ 2m " 
y 1 ' ' ~ 16^ 5 / 2 (mp)3/ 2 ■ 
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FIG. 1: The plot of the ratio F = g ren {x;x) m ^ /g ren (x;x) m =o for v = 0.01 (thin curve) and v = 2 (thick curve). 



In the opposite case, if mp <C 1, the expression for Q ren (x\x) is given by Eq. I|12|) . The plot of the ratio F 
G ren (x;x) m ^o/G ren (x;x) m =o is reproduced in Fig. Qas a function of mp. 
To calculate the energy-momentum tensor we use the following formula |22| 



1 



{Tfx,u) ren = -j Hm Im{Tr[ 7( L At (V„ ) [5 + 5 c ]'- e " - g *\ x ,[S + S*\™ )I(x';x)]}, 



(13) 



4 x' — >x 

where S c is the charge conjugate spinor Green function S c = —C"f^ T S*j^ C" 1 = —"f^S*^^. In above formulas 7^ 
are the gamma matrices in Lorentz regime: 7*- -* = 7j, , 7^ = — vy^ , and r = it. In manifest form 

At* 



S c (x; x) — 



v 



E l A f°° d,S 2 Ar ! +Az ! +p ! +p ,J 



-5 9 c 



2s 3 



2s 



S4) 



(14) 



where 



si; 



5'i 



diag (Ve- <A ^, V , I^e"^, V) , S 2 C = diag (-J^ + e- iA P -I p _ , Ip+e^, I f3 _ ) , 
I ^ + e- ,;A ^ \ 



-i^e"' 7 ^ 






c 

h- 





) 








-e-^{p'l [)+ -pl p _) 
\-e^\p'Ip_- P I 0+ ) 



e-^{p'l 0+ -plp_)\ 
e^'{pfl p _-plp + ) 



(15) 



and the argument of the Bessel functions is z — The manifest form of the spinor Green function S is given 

by Eq. (|SJ). In order to renormalize the spinor Green functions S and S c we have to subtract from them the same 
expression with u = 1. It means that the renormalized Green functions have the same form in which we have to make 
the changes 



vl 



vl, 



P- 



yI P+ 



p+ 



In what follows we assume these changes in the Green functions. Because the renormalized Green functions are finite 
in the coincidence limit we may use 8* instead of g v . To obtain the energy-momentum tensor we use the following 
relations 



5> 



E 



(16a) 
(16b) 
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which can be easy proved by changing Z — > — Z — 1 in the sum in the left hand side. Let us briefly show the main steps 
of calculations (for simplicity we consider the case v < 2). 

FromEq. JTSJ we conclude that the energy-momentum tensor has diagonal form with components (no summation): 

T[! = -- lim ImTr 7 £(^-^)(S + ST e "- 

In this formula the coincidence limit for all variables except has already done. The contributions from S c coincides 
with contribution from S. Let us consider consequently different components of the energy-momentum tensor. 
T t ': There is non-zero contribution only due to S2, because S± and S4 do not depend on the time: 

{\ ren 

T 2 Z : The calculations of the component is similar. The non-zero contribution is only due to S3 and we obtain 
T£ = T}, as should be the case due to boost-invariance along z-axis. 

TP; In the coincidence limit S± and S4 survive. The term Si gives no contribution because 

\un{d p -d pl )e-^ 1 I 0± { P -f) = Q. 
p —>p Zs 

The non-zero contribution appears only due to S4. By using the relations 

1 2 p 2 



and l|l()a|) one has 

lim l P L (d P d p ,)S A = -2I _(!-h P L 2 = 2I p _(!-)h. 
p'—>p zs zs 

Therefore, we obtain = T\. 

T£: In the coincidence limit Si and S4 survive. The component Si has no contribution to T£. Indeed, by using 
Eq. (|16b|) one has 

Tr lim 7 ^(^-^)e iiA ^ 1 =2^_(^-)Tr 7 ^diag(l,-l ) l,-l) = 0. 

tp'—xp Zs 

The non-zero contribution exists only due to S4. One has 

lim {d v - d v ,y lAv S A = 2il lim S 4 + lim {d v -d^)S 4 

As a consequence of the Eq. ()16a|l . the last term gives no contribution to T£ and we arrive at the expression 







<p 

ren 



Thus, using the results of the last section we obtain 



T* = T , I^ = -3T +Ti. 
Therefore one has the following structure of the energy-momentum tensor 

(T„T e " = r diag(l, 1, -3, 1) + T 1 diag(0, 0, 1, 0), (17a) 



where 



,7s ■> f p~ iz — — 8 sin f 
To = -ri-,1 / — - - dz, (17b) 



167T 3 ./n S 



X 



vm z cos ^ f°° Ki (2mp cosh f ) sinh % sinh ^ , 



ir 3 p J Q cosh I cos vy — cos -nv 
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In the massless case, m = 0, we get 



vjh (y 2 -\)(7v 2 + 17) 
8tt 3 p 4 ~ ~ 2880tt 2 p' 
Tx = 



T = = ooon-2„4 » ( 17d ) 



and in the massive one, we have 



^m 2 cos^ [°° K 2 (2mpcQsh.%) sinh | sinh ^ , 



2tt 3 p 2 J cosh 2 | cos vy - cos 7Tf 

The energy-momentum tensor obtained obeys the conservation relation: 



v M r£ = 0. 



Indeed, the non-trivial equation for v — r reads as 



d r T- + -(T r r -T^) = 0. 
r * 

Taking into account the manifest structure of the energy- momentum tensor given by (|17J) . above relations takes the 
following form 

d r T + i(4T -T 1 )=0. 
r 

It is easy to verify that this equation is valid due to well-known relation: 

K' 2 {z)=--K i -K 1 . 

z 

The term T\ appears due to the non-zero mass of the field. It may be represented by the following form 

Ti = -m 2 Tr[G ren {x;x)} = mTi[S ren (x; x)}, 
and gives right the trace of the energy-momentum tensor 

(T[t) ren =T 1 . 

For mp 3> 1, the energy- momentum tensor is exponentially small 

vRt m 4 e- 2m P 



To 



8tt 5 /2 (mp) 5 / 2 



Therefore the energy is localized very close to the string in a radius smaller than the Compton length of the spinor 
particle, p < mT 1 . For mp <C 1 the expression for Tt is given by (|17d|) . The plot of the ratio T , ( m ^o)/^o(m=o) looks 
very similar to the plot in Fig^ 



V. CONCLUSION 



In this paper we considered the Dirac field in the space-time of an infinitely thin and straight cosmic string. We 
found in manifest form the Green function and heat kernel of the squared Dirac operator. We showed that the trace of 
the first heat kernel coefficient is minus two times the corresponding coefficient for the scalar field in accordance with 
Rcf. |20| and found the first three coefficients in manifest form. Thus, the energy- momentum tensor of the massive 
spinor field was obtained in manifest form. In the massless case the tensor is expressed in terms of the second heat 
kernel coefficient. The mass of the field brings an additional parameter, the Compton length of the Dirac particle. 
Due to this parameter the energy-momentum tensor has an Yukawa type dependence over the distance and the energy 
density of the vacuum polarization is concentrated close to the string, in a domain such that p < m _1 . 
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